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PACS 47.63.mf - Low- Reynolds-number motions 

PACS 05 . 65 . +b - Self-organized systems 

PACS 87. lO.-e - General theory and mathematical aspects 

Abstract. - We consider a simple model of an internally driven self-rotating object; a rotor, 
confined to two dimensions by a thin film of low Reynolds number fiuid. We undertake a detailed 
study of the hydrodynamic interactions between a pair of rotors and find that their effect on the 
resulting dynamics is a combination of fast and slow motions. We analyse the slow dynamics 
using an averaging procedure to take account of the fast degrees of freedom. Analytical results 
are compared with numerical simulations. Hydrodynamic interactions mean that while isolated 
rotors do not translate, bringing together a pair of rotors leads to motion of their centres. Two 
rotors spinning in the same sense rotate with an approximately constant angular velocity around 
each other, while two rotors of opposite sense, both translate with the same constant velocity, 
which depends on the separation of the pair. As a result a pair of counter-rotating rotors are a 
promising model for controlled self-propulsion in two dimensions. 
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Soft active systems are composed of interacting units 
t[hat consume energy and generate local motion and me- 
chanical stresses. They show a rich variety of collective 
behaviour, including dynamical order-disorder transitions 
and pattern formation on various scales. A much studied 
example are collections of self-propelled objects [THT5] 
which independently translate due to internally generated 
motions by the objects themselves. 

In this article we study, a particularly simple class 
of self-driven particles and their interactions. These 
are self-rotating objects (that we call rotors) lying in a 
plane [TBUlSj . Our motivation is two- fold. First recent 
developments [T^H^ have shown the possibility of ex- 
perimentally realising such self-driven rotors on the nano- 
scale. In addition recent progress in nano-scale microfab- 
rication techniques have been able to generate synthetic 
chemically driven rotors that rotate fast enough so that 
the effects of their interactions on their dynamics can been 
measured experimentally |23] . Second, a collection of ac- 
tive rotors provide a particularly simple example of an ac- 
tive system as their collective behaviour can be described 
in terms of scalar rather than vector equations. The first 
step in the quest to understand their collective behaviour 
is to develop a detailed understanding of their interactions. 

Motivated by this we investigate analytically and nu- 
merically the dynamics of rotors confined in two dimen- 
sions (the X — y plane) by a viscous film, such as a mem- 



brane. We first introduce the model of a single rotor: 
(which when isolated rotates around its centre but does 
not translate). Each rotor can be described by a scalar a, 
the value of the projection of its angular velocity u> on the 
z direction, that we call spin. We next study the hydrody- 
namic interactions between a pair, in the absence of noise 
when the distance between their centres is large compared 
to their size. We find that the resulting dynamics can 
be described in terms of slow and fast variables that de- 
pend on the relative spin of the pair. We obtain analytic 
expressions for the dynamics of the slowly varying quan- 
tities and find that for opposite spins the hydrodynamic 
interactions lead to a net translation of the pair while for 
like spins they lead to a relative rotation of the pair around 
each other. These leading order analytic results are in 
good agreement with numerical simulations for well sepa- 
rated rotors. Interestingly, our study suggests that a pair 
of counter-rotating rotors could be used to construct a mi- 
cron sized-swimmer in two dimensions, with the speed of 
the swimmer controlled in a precise manner by varying the 
equal and opposite torques acting on the rotors. 

Hydrodynamics of a thin fluid fllm. The ac- 
tive rotors can be confined to two dimensions by being 
adsorbed onto a thin fluid film P^IM^ . The film can 
be described as an infinite incompressible two dimensional 
layer of fluid with (2d) viscosity 77 fllling the plane z = 
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and coupled liydrodynamically to another incompressible 
bulk fluid of (3d) viscosity ije which fills the region z ^ 
(see Fig. [IJ . We consider the fluid dynamics in the vanish- 
ing Reynolds number (Stokes) limit where inertia can be 
neglected [33]. Given a force density F{x,y), the in plane 
equation reads 



F ; V_ 



(1) 



where v[x,y) is the 2-dimcnsional velocity field (in the 
plane 2 = 0), V_L = {dx,dy) is the 2-d gradient operator, 
p{x,y) is the in-plane pressure, := 'r]edzv'\Q± is the 
shear stress of the bulk fluid at the top/bottom of the thin 
film; v'{x,y,z) is the 3-dimcnsional velocity field defined 
in the external region z 7^ which satisfies the Stokes 
equation 



jyeV^v' + Vp' = ; V • = 



(2) 



where V = (Vj^,(?z) is a 3-dimcnsional gradient opera- 
tor and p'{x,y,z) is the bulk pressure. The ratio of the 
two and three dimensional viscosities introduces a length- 
scale I := Tj/'^Ve- Associated with the existence of this 
length-scale are two asymptotic regimes: at lengths r ^ I 
dissipation is mostly due to flow out of plane, and the 
hydrodynamics is similar (but not identical) to that in 
three dimensions; while for r <C I dissipation occurs al- 
most entirely in plane and thus hydrodynamic flow fields 
are quasi- two dimensional. We construct objects out of 
collections of flat disks of radius a subject to point forces 
f at their centres lying entirely in-plane. The viscous drag 
on a disk moving at a constant speed can be obtained 
in an analogous way to the Stokes' solution for a sphere 
in three dimensions and the resulting drag has the form 
[24l[25l[27| 7 = 4:7rr]/g where 17 is a function of I /a. 
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Fig. 1: Thin film. A two dimensional fluid layer characterized 
by a 2d viscosity rj fills the whole plane z = 0. It is surrounded, 
above and below, by another fluid having a 3d viscosity rje with 
which it is coupled hydrodynamically. Objects are constructed 
out of groups of disks of radius a that are subject to forces 
lying entirely in the plane z = 0. 

The interaction between disks embedded in the film can 
be obtained using the Green function H of cq. ([T|), cor- 
responding to the flow, v{r) = H{r — ro) • fo generated 
by an in-plane point-like force, or stokeslet, fo at Tq. This 
description is well suited in the limit where the separation 
between disks is much greater than their radii so they 



can be approximated as point-like objects. The tensor 
/ r cii*^g-»k r (L_^^^ jg ^j^g ^Yiin film equivalent 
and we refer to it as the 



of the Oseen tensor 
hydrodynamic tensor. In the limit r ^ I, 

I r ® r 
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For a set of N discs with forces at their centres, 
(fc G {1, N^) each disk is subject to the dynamic equation 



N 



(4) 



where, for n ^ m Hnm '■= H{r„ — r™), and for n = m we 
have 'Hnn '■= 1/7- This allows us to characterise the flow 
fields on scales much larger than a. 

The rotor model . — A rotor is composed of two 
disks, labelled with the index ?? = 1,2. The disks are 
a fixed distance L apart, oriented along the direction u 
(lying in-plane) . Denoting by x„ the position of the disk 
n, Lu = xi — X2. Equal and opposite point forces (also 
lying in-plane) act at the centre of the disks directed per- 
pendicular to their separation, applying a torque on the 
rotor. The associated force density is fn(5(x — x„), where 
f„ := (— )"+^(T/it-'- so that f := fi = {2- The magnitude 
of the rotor torque, fL and the spin cr = ±1 (the sense 
of rotation), parametrise the rotor. The director is 
obtained from m by a clockwise rotation of an angle 7r/2. 
See Fig. [2] for an illustration. Neglecting thermal noise, 
the equations of motion are xi = (I/7 — H12) • f , ^2 — 
- (I/7 - ^21) • f , where ?ii2 = ?i2i 



27Tri L 




Fig. 2: Instantaneous configuration of a single rotor. For fixed 
/, there are 2 possible choices for the force distribution, de- 
termined by cr, as illustrated, resulting in two different senses 
of rotation. We associate a = 1 with clockwise rotation and 
vice- versa for cr = — 1. 

By construction, the rotor centre R = ^(xi -I- X2) is 
stationary while the orientation axis is not. Its equation 
ii ~ describes rotation with constant angular velocity 
around the centre. To see this it is convenient to introduce 
the normal, z defined as it A = z, a unit vector in 
directed along the negative z— axis. The equation for u 
can then be recast as: 



M = Wo A ti 



(5) 
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with angular velocity, ujq := j^crz. The period of the ro- 
tor is given by Tq '■= 27r/a;o, where wq = 2//(i7). The 
spin a is simply the projection of the angular velocity di- 
rector a>o on to z. Equation ^ determines u{t) up to 
a phase. We choose u, = a; sm{auJot) + y cos{auJot) and 
— xcos{<TUjQt) — y siii{<TUjQt) . 

Dynamics of two interacting rotors. It is 

straightforward to generalize the approach above to deal 
with two rotors, which we label A and B (see Fig. [3]). 
There are now 4 disks with centres at positions, x^, where 
n = 1,2 and a — A, B. It is convenient to describe the 
dynamics in terms of orientations Ua and the centre posi- 
tions Rq of each rotor, defined by Lua '■= x" — X2 ; Rq 

-(x" +X2 ) , respectively. A single rotor cannot propel it- 
self through the fluid. Interestingly, however, net motion 
of the rotor center occurs when another rotor is close by. 
This can be seen in the equations of motion for Rq, and 
Ua that, in absence of noise, are : 
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receives a correction coming from the interaction with the 
other rotor. It is small (of order e^) compared to the lead- 
ing order 0(1) term, which is unperturbed rotation of u 
with frequency loq. It can be shown perturbatively to be 
a time-dependent effect with zero average value over a cy- 
cle to any order of e. Thus, it does not induce a shift in 
Wo [31] • Therefore, 



•ila = XShl{(TaUJot) + y COs((Tc(CJot) + O(e^) 



X cos{craUJot) — y sin(o'aU;oO + O(e^) 



(8) 



Here Hpaqp ■= H(xp — x^) is defined ioi a ^ P and 
r ^ l.a. It is convenient to define 



(7) 



the relative and mean position respectively of the two cen- 
tres of rotors A and B. The magnitude r of the relative 
distance r is the typical separation between the objects. 
For large separations, r 3> it is natural to expand the 
hydrodynamic tensor using the small parameter f.= L/r. 



Fig. 3: Instantaneous configurations of two interacting rotors. 
The two rotors, labelled with the letters A, B, lie in the plane at 
a distance r between their centres. The spin values are (ja = 1 
and as = — 1- 

Angular dynamics. From eqns. it is clear that 
angular and rotational dynamics are coupled. The equa- 
tion for the director Ua shows that the angular velocity 



Since the isolated rotors do not translate, the first non- 
zero contribution to the translational motion is at least 
of order e. Hence to leading order in e, the rotational 
motion affects translational motion of the rotors but not 
vice-versa. In other words, the individual rotor directors 
Ua are fast variables and their positions Rq, are slow vari- 
ables. Consequently, we are now in the position to study 
the effects of the fast degrees of freedom on the (slow) 
dynamics. 

Translational dynamics. The rotor directors Mq, are 
periodic functions of period To. To study the dynamics 
■of the position Rq, on time-scales much greater than the 
period, t ^ Tq, we use a dynamic variant of asymptotic 
homogcnization |32| . This involves performing an asymp- 
totic expansion in powers of the ratio of the short time- 
scale to the long time-scale, (i/To)""'^, while implementing 
the condition that the microscopic dynamics is periodic. 
From this, given an initial configuration: {Rq(0)}, wc ob- 
tain an 'averaged' equation j33| for the velocity over a 
time period To, = dt 9(Rq for the rotor cen- 

tres, which we denote by an over-bar. This is equivalent to 
calculating the displacement in a period Tq, ARq(To) ob- 
tained by integrating cq. ^ over a period To, where orien- 
tations arc given by eq. ([5]) and dividing the displacement 
by the period To. 

As above, it is useful to use the relative position, r 
and mean position, R, as defined in eqn. ([7]). We use 
a subscript for the positions at the instant i = 0, so 
Ro = R(t = 0) and ro = r(t = 0). As one would 
expect the translational dynamics depends on the rela- 
tive spin of the two rotors. To illustrate this we define 
a± := 1/2((T^ ± (Ts): is non-zero for equal rotors and 
zero for opposite rotors while cr_ is zero for equal rotors 
and non-zero for opposite rotors. Then to leading order 
in e wc find: 

(i) for equal rotors {a a = ctb) R is left fixed and r 
rotates around it with a mean angular velocity 



I ■ 

u;i := — / —Ozdt = cjicr+z 
^0 Jo 



(9) 



as shown in Fig. 2] (a). is the angle formed by r with 
the X axis that increases in the counter-clockwise 
sense and wi = [l f e^) / {2TTr]L'^) = {u}^)le^) / {Lg{l / a)), 
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introducing the time-scale Ti := 27r/aji. Here ep = 
L/ra where tq indicates the magnitude of rp. Due 
to the dilute nature of the solution we have Tq <^ 
Ti, hence r rotates much slower than Ua. This is 
in qualitative agreement with previous work [16lll7j 
using different rotor models; 

(ii) for opposite rotors (cr^ = —ub) r remains constant 
on average and R moves orthogonal to r with a 
mean velocity 

— 1 r'^° . 

R — / ndt = -via^i-Q (10) 
^0 Jo 

as shown in Figure U (b) . Here vi = roWi/2 and we 
identify = ^ since r is on average constant. In 
contrast to identical rotors at low Re [16l[T7|, this 
problem has not received much attention. 




Fig. 4: Slow dynamics of the rotor centres, (a) Equal rotors. 
The mean position R is at rest and the dynamics involves the 
relative distance r only. After a period To, we observe a net 
rotation of an angle AS. Here ga = crs = —1 (b) Opposite 
rotors. In this case R moves and r on average stays at its 
initial value. After a period To, there is a net translation of 
R along the direction Vq. Here (Ta = 1. 

We consider the intrinsic reference frame described by 
unit vectors roi '^cf (which are defined in terms of the 
vectors R^i, R^ at t = 0). They satisfy Tq A = z. The 
results above can be expressed in terms of the velocity of 
the centres of each rotor. The slow dynamics of rotor A, 
can be expressed compactly by the velocity 

R^(t) = i(uJrAr(t)) +R, (11) 

with the corresponding expression for rotor B. This rep- 
resents either the effect of a net rotation, occurring only 
when the rotors have equal spins, or a net translation, if 
they have opposite spins. Since for a pair of equal ro- 
tors, r rotates, it is convenient to use the parametrization 

r{t) — To ^ cos r2(t)ro + sin ^^(i)^,^^ where VL{t) :— ujia+t. 

It is interesting to note that symmetry considerations 
alone indicate that the qualitative nature of the motion 
of a pair should be independent of the detailed form of 
rotor [inilTT] . If the separation of the rotors is large com- 
pared to their size, the net flow is well approximated by 



superposing the flow fields generated by the individual ro- 
tors. On timescales long compared to the period of the 
rotation of u, their combined motion can only depend on 
the two directions ro and characterising the configu- 
ration of the rotors. Moreover, the flow field generated 
by each rotor has rotational symmetry. From this it fol- 
lows that any net velocity has no component on the ro 
direction. Hence, it can only have a net component along 
the rQ direction. Implementing this for two identical ro- 
tors, one easily sees that their two centres arc convected 
by equal and opposite flows - leading to rotation; while 
when they have opposite spin they are both subject to the 
same velocity and translate together. This insight can be 
very productively combined with analytic techniques to 
investigate the slow dynamics. 




Fig. 5: Dynamic trajectories of rotor centres in the plane, ob- 
tained from numerical simulations. Continuous lines represent 
the instantaneous positions of the rotor centres incorporating 
both 'slow' and 'fast' motions; solid circles show its slow en- 
velope, at time intervals that are multiples of the period Tq. 
(a) Equal rotors with cta = —1. The trajectories of the rotors 
oscillate in the radial direction while performing a revolution 
around their common centre. The slow dynamics, however, de- 
scribes a circular trajectory, (b) Opposite rotors, with aA ~ I- 
The two rotor centres oscillate around the radial direction while 
moving on average in a direction orthogonal to their separation: 
the 'slow' motion is a straight line in this direction. The ini- 
tial relative orientation is ro = x + y and the parameters are 
I — 5^m, a = 0.5^m, / — lOpN, L — 6^m, ro — 15^m. 

Comparison with numerics. We compare analytical 
results with numerical simulations that involve the full 
tensorial expression of the hydrodynamic tensor, allow- 
ing us to investigate the dynamics of two rotors beyond 
the regime where e 1. Equations © are implemented 
in C and integrated numerically using the Taylor algo- 
rithm [311 . The accuracy of the results can be improved 
by simply using a smaller integration step dt. Here, we 
have used a time-step dt ~ 10~^^. We consider many 
internal cycles To and keep track of positions and angles 
of rotors both i) at each time step (instantaneous) and ii) 
after a period To. We have checked that a different phase 
prescription for u" does not affect the dynamics, which is 
qualitatively unchanged even in the case of a phase differ- 
ence between the directors. 

For e <C 1, results are in quantitative agreement with 
approximate analytical solutions, with which they are 
compared in Figure (O as a function of the initial sep- 
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Fig. 6: Net changes in a cycle as functions of eo- Approx- 
imated analytical results (continuous curves) compared with 
numerical simulations (squares), (a) Equal rotors. The net 
angle \A9\ in a cycle, (b) Opposite rotors. Net displacements 
|A7i| in a cycle, rescaled with the initial separation ro. In both 
cases parameters are the same of figure ((5| (except ro). 



aration tq, at constant L. As e is increased, deviations 
appear as expected, however the qualitative behaviour re- 
mains unchanged. There we plot the net rotation angle 
after a period Tq from simulations and the approximate 
analytic result given by 



Lgil/a) 



(12) 



and similarly the net displacement, rescaled with tq, com- 
pared with the approximate analytic result 



AR 



= TqVi = 



Lg{l/c 



(13) 



Comparison to three dimensional hydrodynamics. It 
is also instructive to consider (briefly) rotors that move in 
an unbounded three dimensional fluid of viscosity 77. Here 
rotors are composed of spheres of radius a, rather than 
disks, linked together. We assume that their directors u" 
can be confined to two dimensions by some mechanism. 
The dynamics of point forces f = fcru^ acting on the 
spheres and lying entirely in that plane can be represented 
by means of the Oseen tensor [3D]: H{v) = ^+^®^ for 
r ^ a. The phenomenology that we find is qualitatively 
similar to rotors moving in a thin film, with only quali- 
tative differences. The centre of a single rotor is at rest 
and its director u rotates with a constant angular veloc- 
ity oj^'^^ = :;^{-^ ^ Tl)- '^^^ interacting rotors show 
again a combination of fast and slow dynamics. The slow 
motion is a rotation around the common centre for rotors 
having the same spin, and a net translation for rotors hav- 
ing opposite spin. These slow motions are described again 
by means of equations that have the same form of eqs (|9]) 
and (|10[). with same directions but different magnitudes. 

To leading order in e these magnitudes are w^'*'' — '^"^ 



for same rotors and 



•jf'^^ = ^rowp''\ for opposite rotors. 
Similar results can also be obtained also by replacing 
the model of rotor with a single spinning sphere of radius 
h with an angular velocity O, of fixed magnitude, directed 
along the ±i axis. The fiow field at distance r, generated 



by a sphere in placed at the origin, is given by 
V = —it A r. 



(14) 



For separations large compared to the sphere radius, the 
velocities of the two sphere centres can be approximated 
by superposing the fields generated by each sphere. Hence 
for rotors of like spin, i.e. flA = ^b, the motion of the 
pair is a rotation about their common centre with angular 
velocity given by A v, with v the same as in eq ([H]). In 
contrast, for opposite spin, each sphere translates in the 
direction perpendicular to the separation vector r, with a 
velocity given by ([T4|) . 

A pair of rotors as a model swimmer. The dynam- 
ical behaviour of two rotors described above suggests a 
natural application. Since two opposite rotors move in a 
straight line, they provide a nice model of a self-propeller, 
whose force density has the form of a torque dipole. The 
mean propulsion speed depends on the rotor separation 
and follows from eqn. (|T3| as ii ~ AR/Tq. The direction 
of motion can be varied easily by controlling the torques: 
the motion can be reversed by simultaneously inverting 
both the torques. The swimmer can be made to change its 
direction by reversing one of them for a short time since 
as discussed above, when the two spins have the same 
value a = 1 { a = —1) the swimmer will rotate clock- 
wise (counter-clockwise) with an angular displacement (in 
1 cycle) given to leading order by eqn. ((T^ . 

Taking account of fluctuations will lead to Brownian 
motion along the direction of the rotor separation with a 
diffusion constant D proportional to the magnitude of the 
fluctuations. This means that the speed of the propulsion 
which will fluctuate around the average -0. As a result, our 
description will remain valid until the separation of the ro- 
tors becomes comparable to rotor size rg ^ L, when the 
rotors strongly interact with each other. This can be used 
to estimate the lifetime of the rotor-swimmer as a first 
passage time off a particle diffusing in one dimension [35] . 
To avoid this, and obtain a permanent swimmer, two ro- 
tors could be mechanically linked together a fixed distance 
apart by a rigid linker [36) . 

Finally we make some general observations about the 
properties of this self-propeller. As required by the scallop 
theorem |37j . a low Reynolds number swimmer must un- 
dergo cyclic shape deformations that are non-reciprocal re- 
quiring at least two degrees of freedom. It is clear that the 
rotations of the pairs of rotors prescribe a non-reciprocal 
cyclic motion. So despite being force-free, and torque-free, 
an assembly made of two rotors can swim, as a result of 
hydrodynamic interactions. It provides another example 
of an idealised simple swimmer model pSMD] that it might 
be possible to construct on the nanoscale. A possible ex- 
perimental realisation would be a pump created by dielec- 
tric spheres made to rotate by a multi-bead optical trap. 
It is noteworthy that more complex rotor-swimmers com- 
posed of two torque dipoles [TB| of opposite sign can also 
be used to construct an albeit less efficient swimmer [3T] . 
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In conclusion, we have investigated analytically and nu- 
merically the dynamics of rotors confined to two dimen- 
sions by a viscous film (such as a membrane) sandwiched 
between a less viscous bulk fiuid. We first introduced the 
model of rotor and subsequently studied the interaction 
between two such rotors. We found that the resulting dy- 
namics depends on the relative spin of the pair and that 
it involves fast and slow quantities. Concentrating on the 
slow variables, we find that for opposite rotors the hy- 
drodynamic interactions lead to a net translation of the 
pair while for like rotors they lead to a relative rotation of 
the pair around each other. There are a number of possi- 
ble experimental realisations of this system. For example, 
qualitative studies of chemically driven 'nano'-rotors have 
already been performed. A promising direction for future 
studies is to confine them in thin films giving greater con- 
trol of their behaviour and allowing a more quantitative 
study of their interactions. Another natural application 
is self-locomotion. Two rotors can be used to construct a 
swimmer: when the torques are equal and opposite, the 
swimmer translates; when the torques are the same, it can 
rotate in a controlled way with respect to its original di- 
rection, as discussed above. Alternatively, two tethered 
rotors with opposite spin can generate a net fiow along 
their common axis, while rotors with same spin would 
generate a fiow with rotational symmetry around them. 
Finally, this study serves as basis and provides a natural 
framework to describe the collective behaviour of many 
rotors, that we will address in the future [?T] . 
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